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Abstract 

o , 

Recently it has been understood that flavor-changing processes mediated by 
Higgs bosons could be a new and powerful tool for discovering supersymmetry. 
In this paper we show that they may also provide an important method for 
constraining the parameters of the minimal supersymmetric standard model 
Qh| (MSSM). Specifically, we show that observation of B s — > at the Tevatron 

O ' implies a significant, model-independent lower bound on tan/3 in the MSSM. 

This is very important because tan/3 enters crucially in predictions and in- 
terpretations of the MSSM, though it is difficult to measure. Within specific 
models, or with other data, the bound becomes significantly stronger. 
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Over the next several years, the Tevatron at Fermilab will be searching for signs of 
physics beyond the Standard Model (SM), by directly producing and observing new 
particles, and by observing rare processes at rates inconsistent with SM predictions. In 
the latter category falls the search for the rare flavor- changing neutral current (FCNC) 
decay B® — > fi + fi~. Observation of this decay at the Tevatron would necessarily imply 
new physics since the predicted rate for this process in the SM is far below the search 
capabilities of the machine. However, supersymmetry (SUSY) naturally predicts large 
enhancements in the decay rate mediated by neutral Higgs bosons [1, 2, 3, 4], and 
in some cases yields branching ratios three orders of magnitude above the SM. So 
observation of the decay would be strong, albeit indirect, evidence in favor of SUSY. 

But observation of B s — > fifi actually tells us more. We will show that it is possible 
to deduce bounds on the fundamental SUSY parameter tan f3, the ratio of the two 
Higgs expectation values, from a signal. This is particularly important because tan f3 
is very difficult to measure — there is no general technique for measuring it at hadron 
colliders — yet almost all SUSY observables depend on it. And a lower limit may be 
almost as useful as an actual measurement because many SUSY observables quickly 
saturate as tan/5 increases. But we can obtain a limit precisely because the rate 
for B s — > /!/! does not saturate. Instead it is strongly dependent on tan/5, rising 
(and falling) as tan 6 (3. As a secondary feature, we will also be able to obtain some 
information on the mass scale of the new SUSY Higgs bosons, since the branching 
ratio scales as the fourth power of their masses. 

We will begin by considering two very general scenarios in which all flavor mixing 
does/does not come from the CKM matrix and in each case we will arrive at a very 
strong bound. But we will also consider some specific models of SUSY breaking and 
show that even more stringent bounds can be obtained for these. Our focus will 
remain on the Tevatron because it has especially good sensitivity to the signal even 
if it does not reach its full luminosity potential. However the B-factories and LHC 
can also use this, and related processes, to further probe tan (3 and the Higgs sector. 
Unfortunately, if a signal is not seen, non-observation cannot be used to draw absolute 
conclusions about the parameter space of the MSSM other than to rule out specific 
model points. It is always possible to choose SUSY parameters in such a way to push 
the signal down to the level of the Standard Model where it would not be observed. 

1 Higgs-Mediated B — ► fifi 

First, a little theoretical background. The question of flavor-changing neutral cur- 
rents (FCNCs) mediated by Higgs bosons was first addressed two decades ago by 
Glashow and Weinberg [5]. It had always been obvious that the Higgs boson of 
the minimal standard model could not have flavor-violating couplings since the cou- 
plings of fermions to the Higgs field defines the fermion mass eigenstates and thus 
also defines our notion of flavor. But in models with two (or more) Higgs fields, 
the fermion mass eigenbasis can be different from the Higgs interaction eigenbasis 
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Figure 1: New contributions to the d-type quark masses in a minimal flavor model [(a)] and 
in general [(a) and (b)]. 



and thus Higgs-mediated FCNCs can occur. In order to avoid large FCNCs incon- 
sistent with experiment, the authors of Ref. [5] proposed several solutions. One of 
those solutions is the "type IF two-Higgs-doublet model: one Higgs field (H u ) cou- 
ples only to up-type quarks, while the other (H d ) couples only to the down-type: 
C = Q l YuUrH u + Q L Y D D R H d . This guarantees the alignment of the Higgs boson 
interaction eigenbasis with the fermion mass eigenbasis. Such a structure can be pro- 
tected from quantum corrections by any number of discrete symmetries under which 
the two Higgs fields transform differently. 

The MSSM possesses the structure of a type-II model classically. But it possesses 
no discrete symmetry which can protect this structure; all such symmetries are broken 
by the /x-term in the superpotential, which must be present to avoid disagreement 
with experiment. And though the type-II structure is also protected by holomorphy 
of the superpotential, this too is ineffective after SUSY is broken. Thus new terms 
are generated in the low-energy Lagrangian of the MSSM of the form Q l YjjUrH\ + 
Q L Y D D R H\. Whether such terms will lead immediately to FCNCs depends on the 
structure of the new Yjj,d Yukawa matrices [1, 2, 3, 4]. 

In Ref. [1] it was shown that there are contributions to Yd which are flavor- 
violating and can have important effects at large tan /3. It is actually quite easy to see 
why. Consider the diagram in Fig. 1(a) in which charged higgsinos propagate inside 
the loop. If we work in a basis in which the down quarks couple diagonally to if J 1 , 
then the up quarks have off-diagonal couplings to proportional to CKM elements. 
In particular, the coupling of to hjjdR is just the bottom Yukawa coupling, y^. But 
the coupling of to s^r is then ytV ts where y t is the top-quark Yukawa coupling. 
Thus the diagram generates a new interaction s^rIII with coefficient proportional to 
ybUtVts- We can rewrite this coefficient as simply y^e, where e includes not only ytV ts 
but also the loop kinematic and suppression factors. (We ignore the phases induced 
by the sparticles in the loops; work including them is under way [6].) 

When H u gets a vev (v u ), it contributes an off-diagonal piece to the fermion mass 
matrix: 

C=(s R b R ) ( m ° ° )( S h L ). (1) 
I \ Vbtv u m b J \b L J w 

This Lagrangian is written in the Higgs interaction eigenbasis of the fermions, which 
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at tree-level is also the mass eigenbasis; we drop the first generation for simplicity. 
The mass matrix can be diagonalized by a biunitary transformation which mixes the 
sl and b^ by an angle sin^ ~ yb€v u /mb. But since = yb(H d ) = UbVa, we have 
sin# ~ etan/3. So although e is one-loop suppressed, the factor of tan/3 can allow 
sm9 to be (9(1). 

At tree-level, a bb — > Higgs transition is possible, with the Higgs decaying to 
leptons, as in Fig. 2. This occurs by exchange of H d . However the coupling of the 
ci-quark sector to H u has shifted the ci-quark interaction eigenstates away from their 
mass eigenstates. In order to replace the interaction eigenstates on the external legs 
with mass eigenstates, we must replace 

b L -> b' L = cos 9b L + sin 9 s L (2) 

which induces a 1>rSl — > /t/t transition through an ii/^ Higgs, Fig. 2. Then the flavor- 
changing amplitude bs — > /x/x is related to the flavor conserving amplitude bb — > /x/x 
by 

^ 6s ^~smM 66 ^. (3) 

There are several noteworthy properties of the Higgs- mediated FCNCs. First, the 
amplitude for B s — > /x/x scales as tan 3 /3 at large tan/3. One factor of tan/5 come from 
sin#, the other two come from the b- and /x- Yukawa couplings which scale as 1/ cos/3. 
Thus the branching ratio scales as tan 6 /3 and provides an incredibly powerful tool 
for constraining tan/3. Second, at large tan/3, the H d Higgs doublet is essentially 
decoupled from the electroweak symmetry breaking. It contains the physical states 
H°, A and iJ ± which all have roughly equal masses. But because Higgs flavor 
changing must disappear in a model with only one Higgs doublet, the FCNC branching 
ratios must decouple as m\. But this does not mean that the effects decouple as the 
SUSY mass scale increases. Rather, in the limit that all supersymmetric masses are 
taken heavy, while A (or H d ) remains light, the rate for Higgs FCNCs approachs a 
finite constant. 

What diagrams contribute to the off-diagonal Higgs couplings? First, even if all 
flavor violation stems from the CKM matrix alone, there is the one-loop higgsino 
diagram of Fig. 1(a) already considered. Models with only this CKM-induced flavor 
violation are known as "minimal flavor violation" (MFV) models. Note that this 
kind of flavor violation has nothing to do with the "SUSY flavor problem." It is 
always present in SUSY because it is generated by the flavor violation in the CKM 
matrix and cannot be eliminated simply by making superpartners heavy or by aligning 
quark/squark mass matrices. 

More generally there are also gluino loop diagrams contributing to the same in- 
teractions (Figure 1(b)), if flavor- violating LL or RR squark mass insertions are 
non-zero. Such an insertion can be generated in minimal models by renormalization 
group running of the squark mass matrices, or may appear in the underlying theory. 
We make no assumptions in our general analysis about the underlying source of this 
flavor- changing in the squark sector, and we refer to this as general flavor violation. 
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Figure 2: Diagram for bb — > Higgs — > /i/i in the interaction (primed) eigenbasis becomes 
bs — > Higgs — > /x/x in the mass eigenbasis, suppressed by sin#, the sz,-&l mixing angle. 

Because this new source of FCNCs emerges from the Higgs sector, it preferentially 
generates processes involving heavy fermions. Thus the channels in which Higgs- 
mediated FCNCs can most easily be observed are B s ^ — > tt and B s ^ — * \x\x. The 
channels with final state muons are suppressed relative to those with final state taus 
by (m M /m T ) 2 but are much cleaner experimentally. The channels with initial state 
B d are suppressed relative to B s by (Vtd/V ts ) 2 . Thus a machine which produces an 
abundance of B s and can cleanly tag and measure muons is ideal for studying this 
physics. This machine is precisely the Tevatron, either in CDF or DO. The current 
limits on Br {B s ^ — > fxfj) are provided by CDF from Run I at 2.0 x 10 -6 and 6.8 x 10~ 7 
respectively. Given enough luminosity, the _B-factories may be able to corroborate 
any Tevatron discovery using Bd — > /i/i or, if a suitable technique is found, Bd — > tt. 

If CDF and DO each record even 2fb _1 of data during Run II of the Tevatron, 
now underway, they can probe the B s — > /i/i branching fraction below the 10~ 7 level. 
(See, e.g., Ref. [8] for a more detailed study of the CDF capabilities.) We will for 
illustration assume that the Tevatron can discover a signal in B — > /i/i if its branching 
ratio is greater than 10~ 7 . Because the branching ratios scale as high powers of the 
input parameters (m\ and tan 6 /3), small changes in the Tevatron capabilities will 
generate infinitessimal changes in our results. And once a signal is discovered, the 
numerical analysis can be redone to obtain more precise limits. 

We demonstrate below that there is a minimum value of tan/5 consistent with a 
signal at the 10~ 7 level, and calculate it. We also consider the interplay between the 
pseudoscalar mass and tan/3, a result which could be important in searches for the 
additional Higgs bosons at future colliders. Both of these analyses will be done in full 
generality within the MSSM. We will not constrain ourselves initially to any particular 
class of MSSM models, such as supergravity- or anomaly-mediated SUSY-breaking. 
This can be done because there are only a limited number of parameters on which the 
branching ratio will depend and these few parameters can be studied without further 
simplifications or assumptions. However in particular classes of models the constraints 
on tan (3 and are stronger and therefore observation of a signal provides even more 
information. We will consider specific SUSY models after doing the general analysis 
and we will find bounds on tan (5 much stronger than in the general case. 
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2 Minimal Flavor Violation: Higgsino Contribution 



Regardless of any details of SUSY-breaking, the higgsino contribution of Fig. 1 (a) to 
Higgs- mediated flavor-changing processes must be present. First we will consider the 
case where only this contribution is present and then generalize in the next section. 
Thus we begin by considering the MFV models. 

The flavor-changing contribution of the higgsino loop diagram is encoded in a new 
dimensionless parameter Kg such that 
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Br (B s - c ^ CB mi s fir Bs ml ro J [ ] [ - J ] . ( 1) 

Here ^qcd — 1-5 is the QCD correction due to running between the SUSY and B s 
scales, and we take f Bs = 220 MeV. Note that the fraction tan 2 /3/ cos 4 /5 approaches 
tan 6 (3 at large tan /3. 

The parameter k~ is calculated to be [1] 



4 v / 27T 2 sin 2 /5^ Jyr ' 



where mj L and mj are the left- and right-handed top squark masses, \i is the super- 
potential Higgs mass parameter, A t is the top-squark trilinear term, and the function 
/ is defined in Ref. [1]. The important thing to know about / is that it is positive 
definite, symmetric in its inputs and f(x,y,z) ~ l/max(a;, y, z) up to a constant of 
0(1). In particular, f(m 2 ,m 2 ,m 2 ) = l/(2m 2 ) and / (m 2 , m 2 , 0) = 1/m 2 . 

In order to maximize the function (and therefore maximize the branching 
fraction one can obtain for a given {m^,tan/9}) we need only consider the four pa- 
rameters m~ , m~ , fj, and A t . This is easiest to do by considering the limits in which 
k~ could become large. 

First consider the limit in which fi is much larger than the squark masses. Then 
~ A t j ji since /i appears both in / and as a prefactor. Thus it appears that Kg can 
grow unabated as A t becomes large. Such a runaway behavior would generate arbi- 
trarily large branching ratios and invalidate our claim to a bound on tan (3. However, 
there is a "cosmological" limit to A t : as A t increases far beyond the squark masses, 
the usual electroweak vacuum become unstable with the true vacuum breaking QED 
and QCD [9]. We apply this constraint by applying the famous condition 

A 2 <3(ml+mf R +m 2 2 ). (6) 

(The last term, m|, is the mass parameter for H u appearing in the Higgs potential. 
Since m| < is necessary to break the electroweak symmetry, we can maximize A t 
by setting m 2 . = 0.) Because the squark masses are much smaller than \i in this limit, 
then A t /(J, can not become much larger than unity at best. A similar argument holds 



5 



100 200 300 400 500 600 700 800 900 1000 

m A [GeV] 



Figure 3: Contours of maximum allowed value for Br (B s — > /i/x) (labelled) as a function of 
{mA,tan P} for minimal flavor violation. For a given branching ratio, the allowed region is 
above the labelled line. 



for large A t but with /j, smaller than the squark masses, except now is more highly 
suppressed, ~ A t 

The actual limit is obtained when \x becomes large along with one of the squark 
masses (with the other as small as possible). Then /(• ■ ■) — > 1/fJ 2 . But the QED/QCD- 
breaking constraint still limits A t , though now it becomes A t < y/3 n so that fiA t f < 
a/3 which implies that 

0.011 
H sm 2 (3 
Plugging into Eq. (4) immediately gives 

„ , „ , 5 x 10- 6 GeV 4 

Br( B ^ W )< <cos6g . (8) 

Thus for itla — 100 GeV (approximately its current lower limit) and a branching ratio 
greater than 10~ 7 , we deduce in an MFV scenario that tan/5 > 11. As uia increases, 
the bound on tan (3 also increases rapidly; for example, if > 200 GeV, one must 
have tan (3 > 18. In general, 



tan (3 > 11 



f rn A \ 3 
UOOGeV/ 



Br (B s -> fjjj) 
1 x 10- 7 



(9) 



where we have taken 1/ cos (3 ~ tan/3 for large tan/3. 

In Figure 3, we plot contours representing the maximal value of the B s — > fifi 
branching ratio consistent with a given choice of and tan (3. Given a measured 
value of the branching ratio, only the region above the line is consistent in an MFV 
scenario. 



6 



3 General Flavor Violation: Gluino Contributions 



In the case of MFV, only the higgsino loop has any significant effect in producing 
flavor violation. However once one moves away from truly minimal flavor violation, 
a new source arises. In a truly general scenario, the g-q-q couplings can be flavor 
changing. In the mass insertion approximation, these violations are moved to the 
squark propagators where they appear as bilinear interactions which mix flavors and 
even chiralities of the squarks. The size of the insertions is parametrized by a dimen- 
sionless quantity, 5^ b = (Am^jj/m^ where a,b are either L or R, and i,j label the 
squark generation. We are only interested in the down-sector 8 2 L \ and 8 R 3 R insertions 
for this work. At present there are no strong experimental bounds on these inser- 
tions and so it is possible that the s~l and b\ (or s R and b R ) are maximimally mixed, 
or 8 2 L \ RR ~ 1. Of course, in particular models a large 8 2 ^ L RR could lead to large 
B — > fT*7 or B — > (f)K s , but constraints from these processes are sensitive to 8 2 \lm 2 ah 
while B — > nil is sensitive to 8 2 % alone. Since we allow very large flavor-mixing, 
we must take care to note that the mass insertion approximation breaks down as 8 
approaches one; we will discuss this issue shortly. 

The new gluino-induced source of flavor changing in general models generates new 
diagrams in B s — > /ifj, with b and s squarks and gluinos in the loop; see Fig. 1(b). 
For our purposes here we can ignore the higgsino contributions since their maximal 
size is smaller than that of the gluino loop if 8f^ LRR ~ 1. We will also assume that 
only one of the two diagrams in Fig. 1(b) dominates the branching ratio; when there 
is data, a fully combined analysis should be done. Specifically we will choose the 8 2 L \ 
diagram to dominate, though our results are identical if the 8 2 RR diagram dominates 
instead. We will discuss more realistic scenarios at the end of the section. 

We parametrize the gluino contribution as 

«5f= ^8? L/ ,M 3 f 2 (ml,ml,m 2 ~ R ,M 2 ) (10) 

where the function f 2 is defined in Ref . [7] . The branching ratio is calculated from k~ 
just as for in Eq. (4). Notice in particular the lack of a Vt. s suppression in k~; this 
will allow much smaller values of tan /3 to be consistent with an experimental signal. 

There are two limits in which we can simplify the above expression. For small 
b L -SL mixing, we can usually take = m~ in Eq. (10) and then work in the mass 
insertion approximation. However, if there is a large hierarchy in between s~l and 
Sl, or if there is large mixing, that approximation breaks down. Instead, it is far 
easier to work directly in the basis of the mass eigenstates. Defining du (i = 1,2) 
as &lcos# + sl sin 9 and its orthogonal combination, we must replace the product 
8f L f 2 {---) inEq. (10) with 

cosflsinfl {f(m 2 ~ Li , ■ ■ •) - /(m^, •••)}< \f(m\ L ^ ■ ■ ■). (11) 
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The inequality follows only from the positivity of the function /. Given more infor- 
mation in the future about the squarks (masses or lower bounds) it may be possible 
to find a lower bound on the second term and thus an improved upper bound on the 
branching fraction. 

Like the MFV case there is a possible runaway behavior in computing the dia- 
grams in Fig. 1(b), allowing for arbitrarily large branching ratios. In the MFV case, 
this arose as A t — > oo, which we cut off by demanding that color-breaking minima 
deeper than the SM minimum not appear. In a general case, this arises as ji — > oo 
for which the color-breaking constraint is useless. However there are equally powerful 
constraints which rely on fine-tuning arguments, which have recently been strength- 
ened [12]. In particular, the /z-parameter appears in the Higgs potential and it is 
the minimization of this potential which must supply the weak scale. There is a 
well-known relation among the /i-term and other Higgs soft mass terms which must 
together generate the scale m z : 



1 2 m H d 
-m 7 = - 

2 z 



m 



tan 2 (3 



tan 2 (3 



(12) 



If \n\ ^> mz then a fine-tuning must be arranged among mu u i rriH d and jj, in order 
to generate m z on the left hand side of the Eq. (12). We will require \n\ < 500 GeV, 
which is a statement that we allow less fine-tuning in the electroweak potential than 
about one part in 60. However in our expressions we will show how to scale our 
results for other choices of /z, in case the reader wishes to apply their own fine-tuning 
constraint. 

In order to calculate the upper bound on the branching ratio, it is useful to note 
that the product M 3 f(M 2 ,m~ L ,m~^) reaches its upper bound when the squarks are 
as light as possible and degenerate, and M 3 ~ 2.1m~. Then one can derive the 
semi-analytic bound (assuming maximal mixing): 



k~ < 0.009 



leading to 



and 



Br (B a -> fifi) < 



500 GeV 
3.0 x 10- 3 GeV 4 



tan (3 > 4 



m A 



100 GeV 



BrfS, 



500 GeV 



tan 6 (3. 



1 x 10- 7 



500 GeV 



(13) 



(14) 
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Contours of the maximum allowed branching ratio are plotted in Fig. 4 as a function 
of tua and tan (3 for the case of maximal mixing. For small mixing we can return to 
the mass insertion approximation in which case 
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(16) 
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Figure 4: Contours of maximum allowed value for Br (B s — > fifj,) (labelled) as a function of 
{rriA, tan 0} for non-minimal flavor violation with 5 2 ^ L = 0.1. For a given branching ratio, 
the allowed region is above the labelled line. 

The bounds presented in this section should be considered the absolute limits 
on tan/? given a signal at the Tevatron at the 1CT 7 level. Nowhere did we inject 
theoretical constraints stemming from some particular model of SUSY. On the other 
hand, if we had a model, or other information, we could devise much stronger bounds. 
For example, many observables in SUSY depend on /i, so that once SUSY is found, 
we will very quickly have constraints on fi, which would allow more precise bounds. 
And if we have some theoretical or experimental reason to believe in a particular 
model of SUSY-breaking, then we know even more. As we will discuss in Section 4, 
models produce bounds which are much stronger than the general case. 

3.1 "Nearly-Minimal" Flavor Violation 

For most purposes, the bounds in the previous section are too general. SUSY models 
which are truly generic and general have an embarassing wealth of flavor-changing 
mass insertions like 8f^ L . Most of these insertions must be miniscule in order to avoid 
dangerously large FCNCs in well-studied processes, such as K°-K mixing. Thus 
one usually imposes on SUSY some form of organizing principle which prevents most 
of these insertions from becoming large, whether that be mass degeneracies or mass 
matrix alignment. 

However even in models which naturally solve the SUSY flavor problem, one 
usually still expects some flavor violation to reappear by way of renormalization 
effects. The most common source of these is the presence of the up-quarks Yukawa 
matrix in the renormalization group equations for the left-handed squarks. Because 
all Yukawa matrices cannot be simultaneously diagonalized, non-zero 5 % l L will be 
generated. (The RR mass insertions will not be generated at leading order.) However 
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there is a significant difference between these mass insertions and those considered in 
the general case above: here 5 l [ L will be proportional to the corresponding element 
of the CKM matrix. For example, 5f? L will be naturally suppressed by V ts . Thus one 
typically finds 8j? L < 0(V ts ) which in turn would force tan/3 > 9 for Br (B s — > fifi) > 
10~ 7 . However this bound is not a true limit - the model can be tuned in order to 
obtain the same branching ratio but with much smaller values of tan (3. 

Thus, while one expects from theoretical arguments that tan (3 must be greater 
than 9 in order to see a signal at the Tevatron, the true mathematical limit remains 
at tan (3 = 4 as found in the previous section, until the origins of SUSY flavor physics 
are better understood. 

3.2 Relation to Other Observables 

The bound on tan (3 in the general case depends strongly on the LL (or RR) mixing 
present in the strange-bottom squark sector. But there are several other key observ- 
ables which also depend on this mixing and which should be measured or constrained 
over the next few years. Key among them are B s -B s mixing, parametrized by AM S , 
and the rate and CP asymmetry in B d — > 4>K S . On the subject of B d — > (j)K s , we 
will not say very much, but refer the reader to Ref. [13] for a full discussion. 

The situation for AM, is more interesting. In the presence of a single, dominating 
LL (or RR) mass insertion, a calculation of AM S yields the approximate formula: 



where (AM S ) SM is roughly 15 to 20 ps -1 . The same formula also holds with the sub- 
stitution LL — > RR. (To obtain this simple result, we have followed the calculations 
of Ref. [15] for B d -B d mixing, appropriately modified, and assumed all squark and 
gluino masses are essentially degenerate at m.) Since 5j? L can be complex, the SUSY 
contributions can either add to or subtract from the SM prediction for AM,. A mea- 
surement of AM S near its SM value would indicate either large fh or small 5j^ L , either 
one of which would suppress B s — > fifi and thus tighten our bound on tan/3. Once 
SUSY partners are discovered, the value of m will be known and a very stringent 
bound on 5 2 ^ L can be extracted from the AM S data, strengthening our bound on 
tan/3 further. 

4 Supersymmetric Models 

In the previous sections, we performed very general analyses in determining what 
could be learned about the SUSY parameter space by measuring Br (B s — > fifi). Now 
we ask the same question, but in the context of specific models of SUSY-breaking, 
namely the three most commonly studied models: minimal supergravity models 




(17) 
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(mSUGRA, also called the constrained MSSM), minimal gauge- mediated (GMSB) 
models, and minimal anomaly-mediated (AMSB) models. Realistic models will nec- 
essarily yield much stronger bounds on and tan (3 for two reasons: first, there is no 
reason for their contributions to FCNCs to be maximal, and second there are often 
correlations with other processes which suppress the flavor-changing contributions. 
Several analyses of B s — * fifi in the context of SUSY models have been completed [10] 
and a more detailed study of this will appear shortly [14]; here we summarize some 
of the important results of this last study. Namely we will find bounds on tan f3 in 
specific models for an observed Br (B s — > /i/i) and discuss their implications. 

For each of the three models to be discussed, we will vary all the free parameters 
of the models and at each "model point" calculate the B s — > \x\x branching frac- 
tion [11]. However we will also apply several phenomenological constraints to avoid 
points which are unphysical. The most important constraints come from the muon 
magnetic moment [16] and the b — > S7 transition [17]. For the former, we will demand 
that the SUSY constributions to (g - 2)„ fall between -5 x 1(T 10 and 57 x 1(T 10 . 
For the latter, we will demand that calculated rate for b — > 57 fall between 2.1 x 1CT 4 
and 4.5 x 10~ 4 . In both cases these are rather wide windows. For the muon magnetic 
moment, a wide window is required by the inconsistency of the e + e~ and r data sets 
which are used to calculate the standard model hadronic contributions to g — 2. For 
b — > 57 we use a wide window because of the uncertainties in the next-to-leading 
order calculation of the SUSY contributions. 

4.1 Minimal supergravity 

Minimal supergravity (mSUGRA) actually contributes fairly efficiently to Higgs- 
mediated FCNCs. It has three key ingredients that help to generate large branching 
ratios. First, mSUGRA models often have large A t trilinear terms, necessary for the 
higgsino contribution. Second, the extended running from the GUT scale down to the 
weak scale generates considerable mixing between the second and third generation 
squarks, allowing for a sizable gluino contribution. Finally, the pseudoscalar Higgs 
does not come out particularly heavy. 

In Fig. 5(a), we show a scatter plot of points in the mSUGRA model space, each 
one representing a distinct choice of tan/3, My 2 , m>o, and A (we take /i > only), 
consistent with all experimental and theoretical constraints. There are several impor- 
tant lessons to draw from the figure. One should immediately note that mSUGRA 
models exist, consistent with all other constraints, that have rates for B s — > \x\x right 
up to the CDF bound; in fact, models were found that passed all other constraints, ex- 
cept the B s — > fifi bound, indicating that B s — > /i/i is already a non-trivial constraint 
on mSUGRA models. 

We also notice that mSUGRA models with a branching ratio greater than 10~ 7 
require tan (3 > 40. Thus observation of B s — > //// at Run II would certainly indicate 
very large tan/5 if mSUGRA is the correct model. 
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Figure 5: Correlation of tan/3 and Br (B s — > fifi) in a representative sample of mSUGRA 
models. 

One also notices another interesting feature in the figure, namely a bifurcation of 
the points at the largest branching fractions. The points in the upper arm all have 
Aq < while all the points with A > lie in the main body of points. This behavior 
is rooted in the b — > 57 constraint. For A > 0, models with very large B s — > /i/i 
also have chargino contributions to b — > 57 which tend to cancel the SM. Thus the 
resulting Br (b — > 37) is too small and we throw the model point out. As the SUSY 
mass scales increase, the cancellation goes away, but the rate for B s — ► /x/x is also 
suppressed. However, the cases with Aq < exhibit a very different behavior. For 
light spectra, these models can have SUSY contributions to b —>■ sj much larger than 
those of the SM, but with opposite sign, so that after adding the two pieces, the 
resulting rate is still roughly consistent with the SM and with experiment. Then as 
the SUSY mass scale is increased, the two pieces start to cancel to zero, ruling out 
the model points. As the SUSY scale is further increased, the SUSY contribution 
decouples and the SM once again dominates. Thus a gap is generated where the 
models predict Br (6 — > 57) ~ 0. 

Thus models with A > can only be observed in Run II if tan /? > 50 while models 
with A < can have tan/3 down to 40. The upshot of this is that a measurement 
of B s — > fifi and an independent measurement of tan/3 could indicate that Aq < 0. 
However, these measurements alone cannot be used to show that A > 0. 

4.2 Anomaly Mediation 

The simplest models with pure anomaly mediation have only one free parameter, 
an overall mass scale, but are not phenomenologically viable. But one can define a 
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Figure 6: Correlation of tan (3 and Br (B s — > /i/x) in a representative sample of AMSB 
models. 

"minimal" anomaly-mediated model which has two free parameters: an overall mass 
scale for the anomaly mediation, and a separate mass scale just for the scalars. We 
considered this minimal model, by varying both of these parameters, as well as tan f3, 
and studying the resulting spectra. In principle, AMSB shares many of the same 
features that make mSUGRA a nice candidate for discovery of B s — > In practice, 
this is only partially true. In AMSB the b — > s'j constraint is much stronger and 
rules out many model points where a large Higgs-mediated FCNC signal would have 
been predicted. We do not find any models with observable branching fractions of 
B s into muons (> 10~ 7 ), consistent with the result of Baek et al. [10]. However, we 
note that this result depends very strongly on the b —>■ s'j calculation; changes in the 
NLO calculation could easily allow for much larger (or smaller) B s — > fifi branching 
fractions. 

A scatter plot of the allowed parameter space is shown in Fig. 6(a). This figure 
again shows a bifurcation into two separate regions, though this time much more 
distinct than in the mSUGRA case. Here there is no fundamental difference between 
the two regions, but the basic explanation is similar to the case of the last section. 
For AMSB models with large tan f3 and intermediate masses, the SUSY contribution 
to b —>■ 57 tends to cancel the SM, predicting a rate which is too small. For large 
masses, the SUSY contribution decouples, but then so does B s — > However, there 
is a region of very light masses where the SUSY contribution is roughly twice the SM 
contribution, but with opposite sign, so that their sum squared is consistent with 
data. This last region is the upper-right region of the figure. There we see that one 
can have Br [B s — > are large as 10~ 7 only if tan (3 > 55. 
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Figure 7: Correlation of tan j3 and Br (B s — ► fifi) in a representative sample of GMSB 
models. 

4.3 Gauge Mediation 

Unlike the previous two cases, one would not expect a large Higgs-mediated FCNC 
signal in models with gauge mediation. Whereas mSUGRA has a large A t and lots 
of running to generate squark mixing, gauge-mediated models have neither. At the 
mediation scale, A t = (so no higgsino contribution) and there is little running unless 
the mediation scale is large (so no gluino contribution). Indeed, we find that it is 
very difficult to generate observable signals at the Tevatron for GMSB models, even 
without invoking constraints from b — > sj. This is best summarized by Fig. 7, where 
no models are found with Br (B s — > /i/i) above 4 x 10~ 8 , making this signal difficult 
if not impossible to observe in Run II. 

To obtain this figure, we varied all the parameters of the GMSB model, including 
the number of messengers, which we allowed to take values of n.5 = 1, 3, and 5. We 
also allowed the messenger scale, M, to go as high as 10 14 GeV. At such large M, 
one would expect to maximize your signal since both A-terms and squark mixing 
benefit from the running between M and rriz- And indeed the maximum signals are 
occurring for large messenger scales. However even with such a large scale, there is 
little signal in this channel compared to mSUGRA or AMSB models. Thus a signal 
at the Tevatron tells us something significant that has nothing to do with tan (3 — it 
would be evidence against gauge mediation itself [10, 14]! 
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5 Conclusions 



Over the next several years, the Tevatron's physics program has a strong opportunity 
to discover new physics in the rare decay B s — > /i/i. And if either CDF or D0 do 
observe this signature, they will have simultaneously provided very strong evidence for 
a supersymmetric world just beyond our current reach. And they will have provided 
an important insight into SUSY-breaking and its mediation mechanism. And they 
will have placed a lower bound on the parameter tan/3, vitally important to future 
studies of SUSY. It is the existence of this lower bound, and its value, that we have 
derived in this work. 

The SUSY parameter tan (3 is one of the most important quantities which needs 
to be measured in a supersymmetric world. In a top-down approach, the calculation 
of tan (3 from first principles would be an important test of some underlying theory. 
In a bottom-up approach, tan (3 is a necessary input to the Yukawa sector and thus 
relevant to everything from Yukawa unification to radiative electroweak symmetry 
breaking. Almost all predictions that can test the form of the underlying theory 
depend on tan (3 in some way, so it must be measured as a first step in studying and 
testing models of SUSY. Yet tan (3 is infamously difficult to measure, particularly at 
hadron colliders or in rare decays — all methods proposed to date are model and/or 
parameter dependent. 

The decay B s — > fifi may improve the situation dramatically because it is unusu- 
ally sensitive to tan/3, providing an opportunity to bound tan/3 if a signal is seen. 
In this paper we showed that some very general theoretical assumptions that do not 
depend on specific models allow us to put significant lower bounds on tan/3 given a 
signal. In the general case with large flavor changing effects allowed (5f^ L ~ 1), the 
bound is given by Eq. (15), leading to tan/3 > 4. This bound increases quickly with 
decreasing flavor violation: for 5f^ L ~ 0.1 the bound is pushed to tan/3 > 7 (Eq. (16)). 
In more typical models in which 5^ L < V ts or V ub) then tan (3 > 9 if a signal is seen 
at the Tevatron. However, the decay B s — > jiji does not require new sources of 
flavor changing to be present; it can be induced simply from the flavor changing al- 
ready present in the standard model CKM matrix. In such stronger bound 
tan/3 > 11 is obtained (Eq. (9)). These bounds are for a B(B S — > /z/z) > 1CT 7 ; we 
have also shown how to scale the bounds as a function of the branching fraction. 
Finally, we find that the bounds are significantly stronger in often studied models 
such as minimal supergravity. 

For many purposes, a lower limit is as good as a measurement since many ob- 
servables (m h , — 2, dark matter signals, etc) saturate as tan/3 increases. There is 
also an upper limit on tan (3 of order 60 from the requirement that the theory remain 
perturbative up to high (unificiation) scales. Thus a lower limit often is tantamount 
to a measurement. 
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